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By neglecting the relative quark momenta in the propagator term, the two-photon and 
two-gluon decay amplitude of heavy quarkonia states can be written as a local heavy quark 
' field operator matrix element which could be obtained from other processes or computed 

with QCD sum rules technique or lattice simulation, as shown in a recent work on r; Ci f, two- 
photon decays. In this talk, I would like to discuss a similar calculation on P-wave \c0 2 and 
Xbo,2 two-photon decays. We show that the effective Lagrangian for the two-photon decays 
O ^ of the P-wave Xco 2 and Xbo 2 is given by the heavy quark energy-momentum tensor local 

operator and its trace, the QQ scalar density. A simple expression for Xco two-photon and 
two-gluon decay rate in terms of the / Xe0 decay constant, similar to that of r\ c is obtained. 

j> • From the existing QCD sum rules value for f Xc0 , we get 5keV for the Xco two-photon width, 

OO ' 

^ ■ somewhat larger than measurement. 

O 

O \ PACS numbers: 13.20.Hd,13.25.Gv,11.10.St,12.39Hg 

O ' 

I. INTRODUCTION 

^ . First of all, I would like to dedicate this talk to the memory of Professor Giuseppe Nardulli, 

who, with great kindness and generosity has initiated the long and fruitful collaboration I have 
with the members of the Physics Department and INFN at the University of Bari. 

In the non-relativisitic bound state calculation [l[ ^| , the two-photon and two-gluon decay rates 
for P-wave quarkonium states depend on the derivative of the spatial wave function at the origin 
which has to be extracted from potential models, unlike the two-photon decay rate of S-wave 
rj c and r\\, quarkonia which can be predicted from the corresponding J/i/j and T leptonic widths 
using heavy quark spin symmetry (HQSS) [3|], there is no similar prediction for the P-wave Xc and 
Xb states and all the existing theoretical values for the decay rates are based on potential model 
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calculations [ll. 144141] . 

Since the matrix element of a heavy quark local operator between the vacumm and P-wave 
quarkonium state is also given, in bound state description, by the derivative of the spatial wave 
function at the origin, one could express the P-wave quarkonium two-photon and two-gluon decay 
amplitudes in terms of the matrix element of a local operator with the appropriate quantum 
number, like the heavy quark QQ scalar density or axial vector current Q"f^ 75Q. We have thus 
an effective Lagrangian for the two-photon and two-gluon decays of P-wave quarkonia in terms of 
heavy quark field operator instead of the traditional bound state description in terms of the wave 
function. This effective Lagrangian can be derived in a simple manner by neglecting the relative 
quark momentum in the heavy quark propagator as in non-relativistic bound state calculation. In 
this talk, I would like to report on a recent work [15[ using the effective Lagrangian approach to 
describe the two-photon and two gluon decays of P-wave heavy quarkonia state, similar to that 
for S- wave quarkonia jij] . This was stimulated by the recent new CLEO measurements of 
the two-photon decay rates of the charmonium P-wave ++ , XcO and 2 ++ Xc2 states. We obtain an 
effective Lagrangian for P-wave quarkonium decays in terms of the heavy quark energy-momentum 
tensor and its trace and that the two-photon and two-gluon decay rates of Xc0,2 and Xbo,2 can be 
expressed in terms of the decay constants f Xc0 and f Xb0 , similar to that for r/ c and 775, which are 
;iven, respectively, by f Vc and f„.. Then a calculation of f Xc0 and f Xb0 by sum rules technique 



or lattice simulation 
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21[ would give us a prediction of the P-wave quarkonia decay 
rates. In fact, as shown below, f Xc0 obtained in [3] implies a value of 5keV for the XcO two-photon 
width, somewhat larger than measurement. In the following I will present only the main results, 
as more details are given in the published paper 




FIG. 1: Diagrams for QQ annihilation to two photons. 



II. EFFECTIVE LAGRANGIAN FOR Xc o,2 -> 77 AND X bo,2 



77 



By neglecting term containing the relative quark momenta q in the quark propagator 



m (Q 2 c , b 



being the heavy quark charge), the P-wave part of the cc — > 77,55 and 66 — > 77,55 amplitudes 
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represented by diagrams in Fig. Q] are 

M(QQ 77 ) = -e Q c>6 ^ _ ^ _ ^ (1) 
with A^v the photon part of the amplitude and the heavy quark part T^ u given by 



A-fo, = -2ei • k 2 e 2fl k lu + 2ei • e 2 k 2fl k lu 

-2e 2 • kiei^k 2u + (fei • ^(ei^j, + e 2M ei„) (2) 
7> = (gi M - (3) 

which can be obtained directly from the following effective Lagrangian for two-photon and two- 
gluon decay of P-wave heavy quarkonia states 

£cs(QQ -»■ 77) = -icxA^Qid p - d M )7„Q (4) 
ci = -e 2 Qc l6 Pi-^2) 2 /4-m 2 Q ]- 2 

With the matrix element of 6q^ lu = Q(d fl — d n)j u Q between the vacuum and Xc0,2 or Xbo,2 given 
by (Q 2 = M 2 ) 

<0\e Q ^\xD> = T M 2 (-g^ v + Q^Q v /M 2 ), 

< 0\0q^\X2 > = ~T 2 M 2 e^. (5) 



The two-photon decay amplitudes are then easily obtained : 

M(X0^77) = -^, t [M2/ ^ ]2 (0) 
**to->TT) = -A& 1g ^ g p (7) 



with T2 = V3Tq from HQSS and 

A = (\)M 2 (M 2 e l ■ e 2 - 2e x ■ k 2 e 2 ■ k x ) (8) 
A 2 = M 2 e^ v [M 2 ei^e 2u - 2(ei • k 2 e 2fM k lv + e 2 • k\(. x ^k 2v 

+ei • e 2 &i Al A;2 i ,)] (9) 

For QCD sum rules calculation or lattice simulation, it is simpler to compute the trace of the 
energy-momentum tensor Qq^ given by 2mQQQ. We have then 

v{P2)T^u(pi) = 2mQ v(p 2 )u(pi) (10) 
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The problem of computing the two-photon or two-gluon decay amplitude of Xc0,2 and Xbo,2 states 
is reduced to computing the decays constants f Xc0 and f Xb0 defined as 

< 0\QQ\ X o >= rn Xo f Xo (11) 

Thus To is given directly in terms of f Xo without using the derivative of the P-wave spatial wave 
function at the origin. 

T = ?f (12) 

Thus by comparing the expression for XcO and r\ c we could already have some estimate for the 
XcO two-photon and two-gluon decay rates. For / Xc0 of 0(f Vc ), one would expect r 77 (x c o) to be in 
the range of a few keV. 

The decay rates of Xc0,2, Xbo,2 states can now be obtained in terms of the decay constant f xo . 



We have : 



r 77 (xc2j 



{M Xc0 +bf 
( 4\ AirQtal m Ml fl 
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(M Xc2 + by 



(13) 
(14) 



where -Bo = 7r 2 /3 — 28/9 and B 2 = —16/3 are NLO QCD radiative corrections 
This expression is similar to that for rj c : 
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24] 



47rQia 2 em M v Jl 



T 77 (rj c 



The two-gluon decay rates are : 
'2\ 47ra 2 M 3 f 2 



a s (20 - vr 2 



7T 



c0) 



T gg(Xc2) 



V (M Xc0 + by 

4\/2\4vra 2 M3 c2 / 2 



[l + C (a B /n)], 



[l + C 2 (a s /7r) 



(15) 



(16) 



(17) 



15J\9J {M Xc2 +bY 

where Cq = 8.77 and C 2 = —4.827 are NLO QCD radiative corrections. As with the two-photon 
decay rates, the expressions for two-gluon decay rates are similar to that for r\ c : 
'2\ 47ra 2 M r? J 2 c 



r gg(Vcj 



9 J {M Vc +bf 



1 + 4.? 



.a. 



IT 



(18) 



In a bound state calculation, using the relativistic spin projection operator 
are given by 



25, 



M , f Vc and /- 



frjc ' 

fxo 



T K (0) (4 mo) , 

32 7r m% 



12, 



(8vrmQ) 



KB 

M 



(19) 

(20) 
(21) 



which gives the decay amplitudes in agreement with the original calculation 

Comparing with /„ c , we have 
^i(O) 



/Xc0 = 6 



(22) 



which becomes comparable to f Vc . 

Thus by comparing the expression for XcO and r\ c we could already have some estimate for the 
XcO two-photon and two-gluon decay rates. For / Xc0 of 0(f Vc ), one would expect L 77 (x c o) to be in 



the range of a few keV. As s 
rules value f Xc0 = 357 MeV 



rown in Table 1, the predicted two-photon width of XcO from the sum 



181 ] is however almost twice the CLEO value, but possibly with large 



27] 



theoretical uncertainties in sum rules calculation as to be expected, while a recent calculation 
implies a larger decay rates for XcO- The measured ratio r 77 (x C 2)/r 77 (xco) is then « 0.24 ± 0.09, 
somewhat bigger than the predicted value of about 0.14 as shown in Table 1 together with the CLEO 
measurement of the decay rates [n| which gives (2.53±0.37±0.26) keV and (0.60±0.06±0.06) keV 
for XcO and Xci respectively. 



Reference 


r 77 ( Xc0 )(keV) 


r 77 (x c2 )(kev) 


p r 7T (xc2) 
n ~ rwxco) 


Barbierifl] 


3.5 


0.93 


0.27 


Godfrey[4] 


1.29 


0.46 


0.36 


Barnes [5] 


1.56 


0.56 


0.36 


Gupta[7] 


6.38 


0.57 


0.09 


Miinz[8] 


1.39 ±0.16 


0.44 ±0.14 


u - oz -0.12 


Huang [9] 


3.72 ± 1.10 


0.49 ±0.16 


o n+ 011 

u - lo -0.06 


Ebert[10] 


2.90 


0.50 


0.17 


Schulerfll] 


2.50 


0.28 


0.11 


Crater [12] 


3.34-3.96 


0.43 - 0.74 


0.13-0.19 


Wang [13] 


3.78 






Lavertvfl4l 


1.99-2.10 


0.30-0.73 


0.14-0.37 


This work 


5.00 


0.70 


0.14 


Exp(CLEO)fl6] 


2.53 ±0.37 ±0.26 


0.60 ±0.06 ±0.06 


0.24 ±0.04 ±0.03 


Exp (Average) [161 


2.31 ±0.10 ±0.12 


0.51 ± 0.02 ± 0.02 


0.20 ±0.01 ±0.02 



TABLE I: Potential model predictions for Xc0,2 two-photon widths compared with this work. 



The two-photon Xc0,2,x' C 2 branching ratios are independent of /- 



Xci) 



^(Xc0,X. 



c'O 
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77) = g Q 
6 



4 "em 



a. 



1 + (B Q - Co) — 



7T 



B(Xc2, XC2 -> 77) = 7 Qt (l + (B 2 - C 2 A 
5 at V 7T 



(23) 
(24) 
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with B = tt 2 /3-28/9, B 2 = -16/3, C = 8.77, C 2 = -4.827. Apart from QCD radiative correction 
factors, the expressions for branching ratios are very similar to that for r\ c and r( c : 



B(Vc, ric> -»• 77) = i Qt^ 1 - 8-2 - (25) 

with a s evaluated at the appropriate scale. 

For a s = 0.26, B(r) c — >■ 77) = 3.6 x 10 -4 to be compared with the measured value of (2.8 ±0.9) x 
10~ 4 jlTjj , but this prediction is rather sensitive to a s , for example, with a s = 0.28, one would 
get B(j] c — > 77) = 2.95 x 10 -4 , in better agreement with the measured value of (2.41q'q) x 10 -4 
and for Xc0,2, the predicted two-photon branching ratios would be 3.45 x 10~ 4 and 4.45 x 10 -4 
compared with the measured values of (2.35±0.23) x 10 and (2.43±0.18) x 10 -4 , for XcO and Xc2 
respectively. The predicted branching ratio for Xc2 is rather large and one would need a s = 0.36 
to bring the predicted value closer to experiment. 

Recently, the Z(3930) state above DD threshold found by Belle 28] with mass (3928 ± 5 ± 
2) MeV and width (29 ± 10(stat) ± 2(sys)) MeV, consistent with x' C 2i seems to be confirmed by 
the observation of a similar state by BaBar [29J], with mass (3926.7 ± 2.7 ± 1.1) MeV and width 



(21.3±6.8±3.6)MeV. Belle 



28] gives T 77 (xc 2 ) x B{DD) = (0.18±0.05±0.03) keV while Babar [2 



gives r 77 (x' c2 ) x B{DD) = (0.24±0.05±0.04) keV for this state. If taken to be the 2P excited state 
X' c2 and assuming B{DD) « 0.70-1 [3o||32|, one would get r 77 (x^) = (0.18-0. 24±0.05±0.03) keV 
. This implies ~ (195 - 225) MeV and r gg (x' c0 ) in the range (5 - 10) MeV. 

F° r Xbo,2 potential model calculations similar to that for Xc0,2, gives the two-photon width 



about 1/10 of that for % , which implies f Xb0 = /n 6 /3, smaller than Cornell potential 
/x w =0.46/, 6 . 



331 ] value 



III. REMARK ON THE rf c TWO-PHOTON DECAYS 



Since the predicted two-photon branching ratios for XcO,2,x' c o,2 and for rj c ,7]' c are similar and 
independent of the decay constants, apart from QCD radiative corrections, as seen in Eq. (|23M24p 
and Eq. (|25p . one expects a large two-photon decay rates for rj' c , it would be relevant here to mention 
the problem of the r]' c — > 77 decay rate [3, 35]. The small value of r 77 (r/^) = (1.3 ± 0.6) keV given 



previously by CLEO 



34| is obtained by assuming B(r)' c — > KsKtt) rs B(rj c — > KsKtt). However, 



with the recent BaBar measurement of the ratio [36] 

B{B+ -»• r/ c (2S)K + ) x B(rj c (2S) -> KKir) 



R( Vc (2S)K+ / Vc K + ) 



B(B+ r] c K+) x B(r] c KKir) 
0.096lg;g^(stat) ± 0.025(syst) 



(26) 
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and the Belle measurement |3 



B(B+ -> rj c K + ) x B{rj c -> KKtt) = (6.88 ± 0.77±g;||) x 10 



(27) 



BABAR obtains [M] 



e(??c -> K s Ktt) = (1.9 ± 0.4(stat) ± l.l(syst))%. (28) 

MS quoted by CLEO Q. This new BABAR value for B(rL — > KsKtt) is considerably smaller than 
the corresponding value B(r\ c — > KsKtt) = (7.0 ± 1.2)% [17] for r} c . 

Thus with the BaBar result for B(r/' c — > KsKtt) and the CLEO measurement 

T 77 (^) x B(rf c -> KsKtt) 



one would get 



r 77 (r? c ) x fi(r? c -»■ 2T5AV) 



0.18 ±0.05 ±0.02 



(29) 



r(r? c ' -> 77) = (4.8 ± 3.7) keV 



(30) 



in agreement with the predicted value 



r(r? c / -»• 77) = (4.1 ± 2.3) keV 



(31) 



while the assumption of near equality of the KsKtt branching ratios for rj c and rj' c 



B(rf c -> KsKtt) » B(r/ C KsKtt) 



and the Belle ratio 39J 

B(S -»■ Krj c (2S) x B(r) c (2S) -»• KsK'tt^ 



R(rf c K/ Ve K) 



B(B° -> K?? c ) x e(r? c ->• K s K-tt^ 



0.38 ±0.12 ±0.05 



would lead to 



38] 



(32) 



(33) 



r 77 (r^) = (1.3±0.6)keV 



(34) 



which is rather small compared with the predicted value given in Eq. (|3ip above. 
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IV. CONCLUSION 

In conclusion, we have derived an effective Lagrangian for Xc0,2 and xto,2 two-photon and two- 
gluon in terms of the decay constants f Xc 60 , similar to that for rj c ^ in terms of f Vc b . 

Existing sum rules calculation, however produces a two-photon width about 5 keV, somewhat 
bigger than the CLEO measured value. It remains to be seen whether a better determination 
of f Xc0 from lattice simulation or QCD sum rules calculation could bring the Xc0,2 two-photon 
decay rates closer to experiments or higher order QCD radiative corrections and large relativistic 
corrections are needed to explain the data. 

The problem of two-photon width of rj' c would go away if more data could confirm the small 
BaBar value for B(r]' c — > KsKir) compared with B(j] c — > KsKir). 

As relativistic corrections should be small for P-wave bottomia Xbo,2 states, two-photon and 
two-gluon decays could provide a test of QCD and a determination of a s at the the mi, mass scale. 
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